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There exist methods to reformulate in an exact way the many-body problem of interacting bosons 
in terms of the stochastic evolution of single particle wave functions. For one such reformulation, the 
so-called simple Fock scheme, we present an elementary derivation, much simpler than the original 
one. Furthermore, we show that two other schemes, based on coherent states of the matter field 
rather than on Fock states, lead to an infinite statistical uncertainty in the continuous time limit. 
The simple Fock scheme is therefore, up to now, the only one that was proved to lead to a convergent 
Monte Carlo simulation scheme at all times. 

PACS numbers: 05.30.Jp, 03.75.Fi, 02.70.Ss 

INTRODUCTION 

Consider a gas of N indistinguishable bosonic spinless particles with binary interactions, a situation typically 
encountered in Bose condensed atomic gases 0]. The time evolution of the many-body state vector of the gas is 
described by the Schrodinger equation 

= (1) 

at in 

where Ti. is the Hamiltonian of the system, including kinetic energy, trapping potential energy and interaction terms. 
In the general case of a large number of particles and a large number of modes of the atomic matter field, the Hilbert 
space of the system is far too large for a direct integration of (|l|) to be feasible on a computer. A possible way out, 
intensively explored in particular in the quantum optics community, is to reformulate the evolution equation in terms 
of the stochastic evolution of classical fields, which are much smaller objects. Historically, the first suggestion for 
such a reformulation is based on the Positive-P distribution |j . Other schemes have been proposed more recently, 
based on different representations of the many-body density operator ^, |^ . Generalizations to fermionic systems 
are also possible [Q, ^. 

As these methods are stochastic in nature, their prediction for a finite number of realizations deviates from the exact 
result by a fluctuating, zero-mean quantity, called the statistical error. One of the challenges of such reformulations 
is to prove that the average of the absolute value of the statistical error remains finite at all times, so that an actual 
Monte Carlo simulation converges to the exact result in the limit of an infinite number of realizations ^. To our 
knowledge, a proof of such a convergence property for the multi-mode Hamiltonian problem defined above was given 
only for the so-called simple Fock scheme in [ ij and for its fermionic counterpart |8| . 

The goal of the present paper is twofold. First, it presents a much simpler derivation of the simple Fock scheme than 
the original one in Q: the main idea here is to work directly with state vectors It/i) rather than with the many-body 
density operator. Second, it addresses the issue of the convergence of two other schemes based on coherent states of 
the bosonic field: the simple coherent scheme of Q and a new scheme that we call the Bargmann scheme. It is found 
that these two schemes based on coherent states lead to an infinite mean of the absolute value of the statistical error 
for any finite evolution time. 



^ Let us call 5 the absolute value of the statistical error (to be defined more precisely later on). Let us assume that the probability 
distribution for 5 has a long tail scaling as 1/(5^+^ for large &. The most favorable case is /i > 2, in which case 5 has a finite variance and 
the statistical error for a sample of A/V Monte Carlo realizations scales as Xj yJNr- In the intermediate 1 < ^ < 2 case, the convergence 
is slower, scaling as l/A/'r ^^'^ . In the worst < 1 case, the mean value of 5 is infinite and the statistical error of the Monte Carlo 
simulation grows to infinity with the number of realizations as Mr "^^^^ ^ . For more details, see |p|. 
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THE MODEL HAMILTONIAN 



A model for an ultracold trapped interacting Bose gas in D dimensions with short range interactions can be obtained 
in a second-quantization formalism by using the Hamiltonian 

n = ^AF^'t(r)/io*(r) + AV^4'^(r)«'^(r)*(r)*(r); (2) 

r r 

the spatial coordinate r runs on a discrete orthogonal lattice of M points with periodic boundary conditions; V is the 

2 

total volume of the quantization box and AV = V/M is the volume of the unit cell of the lattice, ho = + is the 
one-body Hamiltonian in the external trapping potential C/cxt (r) , m is the atomic mass and interactions are modeled by 
a two-body discrete delta potential with a coupling constant go. The field operators ^(r) satisfy the Bose commutation 
relations [*I'(r), ^t(r')] = S^.t'/AV and can be expanded on plane waves according to ^'(r) = X^k '^kC*'^'"/ a/V^ with k 
restricted to the first Brillouin zone of the reciprocal lattice; in the plane wave basis, the kinetic energy term has the 
diagonal form %^aj^ak. In order for the discrete model to correctly reproduce the underlying continuous field 
theory, the grid spacing must be smaller than macroscopic length scales like the thermal de Broglie wavelength and 
the healing length. 

THE SIMPLE FOCK SCHEME 

Consider the N particle bosonic Hartree-Fock state defined as 

|iV : 0) = ^a;^|0) (3) 

where the operator aj, creates a particle in the non-necessarily normalized single-particle wave function (j){r): 

a; = 5]Al/0(r)r(r). (4) 

r 

The challenge is here to find an equation of motion for 0(r) in order for the Hartree-Fock ansatz to give an exact 
solution to the full many-body dynamics, i.e. we aim to obtain a dynamics for (f>(r) such that 

d\N : ^) = ^n\N : 0) (5) 
in 

in a sense to be defined. In the following we proceed to evaluate both sides of (H). 

The action of the Hamiltonian (H) on |7V : 0) can be worked out by splitting the field operator in its longitudinal 
and orthogonal components with respect to (fi: 

*(r) = ^a^ + *^(r) (6) 



where the squared norm of (/) is 



(7) 



with the following scalar product 

(</>|^)=^AF0*(r)V(r). (8) 

r 

From this splitting of the field operator it follows that 

^(r)|7V : 0) = \/iV0(r)|iV -1:0) (9) 
■^Hr)\N : (f>) = VWTT^^-^\N + +^{{r)\N :(!>). (10) 
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Inserting these expressions into the Schrodinger equation for the Hartree-Fock state and collecting the terms containing 
the different powers of we finally obtain 



n\N ■.(j)) = VN (^AVx{r)i'Hr)^ \N - 1 : 



^^''^ (^Ayffo</'(r)0(r)vI'^(r)^F^(r)J |7V-2:0) (11) 



with 



X(r) = 



, 9o{N-l) ,^,„„2 go(iV-l) 
"0 H o — 0(r) 7— 

Uf 2 ll^f 



(r). 



(12) 



On the other hand, the variation of a Hartree-Fock state following a variation of 0(r) can be calculated by replacing 
(j){v) by 0(r) + d(p{r) in (^) and expanding with the binomial formula. Up to second order in the variation (i0(r), this 
leads to 



\N : (t) + d(l)) = \N : (j)) + 



^^ ^^ ^) I J2 d(j3{Y) d(l){r') ¥{r) *t(r') | - 2 

r,r' 



(13) 



In the usual Hartree-Fock mean-field approximation, has a deterministic evolution with d(j) cx dt, so that the last 
term on the right-hand side of (^3|) is negligible in the limit dt ^ and therefore the last term on the right-hand side 
of (p^) cannot be accounted for. The best that one can do in the mean-field approximation is to set 



(14) 



thus recovering a nonlinear Schrodinger equation for cj) of the same form as the usual Gross-Pitaevskii equation 
The main idea to recover the exact evolution is to include in d(f> a term dB proportional to Vdt: 



dm = Sx(r) + dB{v) 
in 



(15) 



so that the quadratic term in d(f) indl3) is no longer negligible. The price to pay for this is the apparition of an 
extra term proportional to \/dt in (|13D which has no counterpart in (|l^). The trick is to make dB stochastic with a 
vanishing mean value. Then the average variation of the Hartree-Fock ansatz exactly matches the one required by 
the many-body Schrodinger equation: 



dt 



d\N : (j)) = —n\N : 0) 
in 

provided that the noise term has the following correlation function 



dt 



dB{v)dB{Y') = -QrQr 
in 



|^5r,r'</'(r)</.(r') 



(16) 



(17) 



where the projector Qr = 1— |0)(0|/||()()||^ projects orthogonally to (j). The zero mean property of the noise makes the 
terms proportional to ^fdi vanish in the averaging over noise. The terms of order three or more in dip in the expansion 
( p^ ) are negligible as compared to dt in the limit dt — + 0. In mathematical terms, this procedure corresponds to 
assuming that the evolution of (j) is governed by an Ito stochastic differential equation jl^. Figj^ gives a geometrical 
interpretation of the Monte Carlo sampling of the many-body state vector 1-0) by Hartree-Fock states. The exact 
many-body state vector at any time i > |?/'(i)) is equal to the statistical average of Hartree-Fock ansatz: 



\m)^{\N:m)) 



(18) 
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FIG. 1: Schematic view of the Monte Carlo samphng of the exact many-body state vector |?/)) at time t (cross) by a coUection 
of random Hartree-Fock states : ^h) (stars). In the hmit of an infinite number of realizations, the center of mass of the 
distribution of the stars exactly coincides with the cross. In the TV-body Hilbert space, u is an arbitrary unit vector and c„ is 
the component along u of and of the |A'' : 0^). 



provided that this is the case at the initial time t — Q. 

In practice, a noise with the correlation function ( p^ can be obtained as 



k>0 



■i(k-r+0(k)) 



+ C.C.) 



(19) 



this expression was already discussed in and is easily generalized to finite range interaction potentials. The index 
k is restricted to a half space and to the first Brillouin zone. The ^(k)'s are random angles uniformly distributed in 
[0, 27rl. An alternative choice can be 



dB(v) = 



ihAV 



Qr [0(r)rfe(r)] 



with the d^{r) independent zero-mean random variables with rf^(r) d^{r') = 5r,r' dt. In both cases 

9odt ,|,,|2 



\\dBf < 

We can then calculate the evolution of the squared norm 

d 



^ = WdBlf < 



hAV ' 

2 

and check that it remains finite at all times: 
90 dt 2 



hAV 



the absence of terms linear in dB being due to the orthogonality of dB with respect to (f>. This leads to 

M\t)<\\<l>\\\0)e^, 



(20) 



(21) 



(22) 



(23) 



so that (j) remains finite at all times. The present stochastic wave function reformulation of the bosonic many-body 
problem is exactly the same as the so-called simple scheme with Fock states introduced in [|| from a different and 
more general point of view. Among the schemes based on a Hartree-Fock ansatz for the many-body density matrix 
of the form 



(24) 
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with (j)i,2 evolving according to stochastic differential equations, the simple Fock scheme was shown to be the one 
which minimizes the growth rate of the statistical variance of a around the exact many-body density matrix p(t): 

VarH ^ - ait)\f) = (Xr [cTHt)ait)] ) - Tr [pitf] , (25) 

where ( ■ • • ) represents the average over all stochastic realizations. The last term in the above equation, the so- 
called purity of p, is constant for Hamiltonian evolutions. If the stochastic wave functions and (f>2 are statistically 
independent at the initial time t = 0, they remain so at any later time since their evolutions are independent. The 
mean squared norm of a is then 

(Tr[at(t)a(t)])^(A)2(i) (26) 

where 

A{t)^{N:mN:m^m''' ■ (27) 

A geometrical interpretation of A can be put forward in the spirit of figj^: the mean value of A is related to the mean 
value of the squared distance in Hilbert space between the Monte Carlo ansatz and the exact many-body state vector 



(A) - 111^)11^ -(|||A^:0)-|V^)| 

For the present Fock scheme, it follows from (|2^) that 

{A){t) < {A){0)e^3ot/hAv. 



(28) 



(29) 



this upper bound was already discussed in [Q. In the model of this paper, the number of modes of the bosonic field is 
finite, so that the inequality (29) guarantees that the simulation can predict the expectation value of any observable 
with a finite statistical dispersion. The Monte Carlo statistical variance AO^ on the quantum expectation value 
Tr[pO] of any observable O is indeed limited from above by 



AO^ = (|Tr[crO] |2) - Tr[pO]2 < Tt[0^] (Tr[(7V]) - Tr[pO] 



(30) 



where the traces are taken in the A^-particle subspace and the trace Tr[0^] is a finite quantity for a finite number of 
modes. 

This formalism can be extended to the imaginary-time evolution which allows one to sample in an exact way 
the thermal equilibrium state at a given temperature T: 



This requires solving the imaginary-time evolution equation 



(31) 



(32) 



for a "time" r going from to /? = l/ksT. The use of the projector onto the N particle subspace Vn as the initial 
condition allows one to calculate the physical quantities in the canonical ensemble, i.e. for a fixed number N of 
particles. As we shall prove in full detail in the next section, the projector operator Vn has a simple expression in 
terms of Hartree-Fock states. 

The correct imaginary-time evolution of the Hartree-Fock state \N : <j>) is recovered if the stochastic wave function 
(j){r) evolves according to 



#(r)=-y 



with a noise term dB{r) equal to 



9o{N- 



9oiN- 



{r)+dB{r) 



dB{v) 



drgQ 



i(k-r+e(k)) 



+ C.C.) 



loo 



(33) 



(34) 



We have recently applied this method to study the probability distribution of the number of condensate particles in 
a ID trapped Bose gas for temperatures both above and below the transition temperature 
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PROOF OF (OVER)-COMPLETENESS OF THE HARTREE-FOCK STATES 

In the previous section we have shown how the iV-boson time-dependent problem can be exactly reformulated in 
terms of the stochastic evolution of single particle wave functions; in particular, the thermal equilibrium state in the 
canonical ensemble can be obtained by means of an imaginary-time evolution with the projector Vn as the initial 
state. Here, we shall prove that the projector Vn has a simple expression in terms of Hartree-Fock states and that any 
state vector can be expanded onto a set of Hartree-Fock states. This ensures that the identity ( |T^ ) can be satisfied 
at time t = for any initial many-body state vector |'0(O))- 

It is a well-known fact of quantum optics that the Glauber coherent states form an overcomplete basis of the 
Hilbert space of an harmonic oscillator. In the general multi-mode case, this property can be expressed as 



1 



I?</)|Glaub : 0)(Glaub : (/.| (35) 



where the Glauber coherent state is defined as usual as 

IGlaub:^/.) =e-ll^ll'/2e5:.AV0(r)*t(r)|Q^^ ^3^^ 

and 1 is the identity on the bosonic subspace of totally symmetric wave functions. The integration in ( |35| ) is performed 
over the 2M dimensional space of complex-valued functions defined on the Af-point lattice: 

= Jl dRe0(r) dlm(/.(r). (37) 

r 

Using the well-known formula 

T'jvlGlaub : 0) = -L.e-ll*ll'/2|^ . (38) 
to project onto the N particle subspace, an explicit expression for Vn in terms of Hartree-Fock states can be obtained 

Vn = VnIVn = ^ j 25(/)pAr|Glaub : (^)( Glaub : (j)\VN = C jvlj) \N : lj)){N : (39) 



The constant C is equal to 



t^-^/xiwi iwi"'""-^-"*"'- "^:"^,'" (40) 



The normalized wave function 0(r) is defined as 0(r) = 0(r)/ |10|| and the integration over 4> is performed on the unit 
sphere \\(t>\\'^ = 1 according to the measure Vcf) defined as 

v^^Uf""'-' dU\\v4>. (41) 

As discussed in [ ]TT| , the expression ( |3l| ) is the starting point for an exact numerical calculation of the thermodynamical 
equilibrium properties of the interacting Bose gas. 

Completeness of the set of Hartree-Fock states is an immediate consequence of the result (|3^): any state vector 
belonging to the N particle subspace can in fact be expanded over Hartree-Fock states according to 

\iIj)=Vn\'^)=C jv^{N ■.4>\'^)\N -.4,). (42) 

As required for practical implementation of the Fock scheme, the weight factor can be made real and positive for all 
the (f) just by reabsorbing the phase of {N : into that of the wave function 0: 

{N:^\^) |7V:r(0)(^) (43) 



1^) = C / 

with the phase factor r{(f)) defined according to 



K0)- = S^l^ML, (44) 
{N : (AlV) 
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Since the scalar product of two Hartree-Fock states is equal to 

(iV:02|A^:0i) = (02|0i)'^, (45) 
this family of states is not an orthonormal basis, but rather an overcomplete set. 

STOCHASTIC SCHEMES WITH COHERENT STATES: AN INFINITE STATISTICAL UNCERTAINTY 

The Bargmann scheme 

The derivation of the Fock scheme can be reproduced for the case of unnormalized coherent states, the so-called 
Bargmann states: 

|Barg:0) = Hq e^.^^-^W^^Wjo), (46) 

where cj) is a. stochastic dynamical variable and Hq is a time-independent amplitude. Thanks to (^), one can show 
that any initial state vector \ip{0)) can be obtained as a statistical average of Bargmann states of the form (^); the 
initial value of both (jj and Ho can then vary randomly from one Monte Carlo realization to the other. 

The exact time-evolution given by the many-body Hamiltonian (||) can be shown to be exactly recovered if the field 
amplitude (j){r) evolves according to 

d0(r) = 5/io0(r) + dB{r) (47) 
in 



with a noise correlation function equal to 



dt go 



dB{r) dB{r') = -r ^ 0(r)0(r'). (48) 
In practice, we take the following implementation of the noise which minimizes the squared norm of dB B: 



dBir) = J||£<^(r)^(eK'^-+^e^)) -fee), (49) 

* k>0 

where k is restricted to a half space and to the first Brillouin zone. Within the framework of [Q, this Bargmann scheme 
recovers one of the possible solutions of the consistency equations in the case of a coherent state ansatz, namely the 
one with a constant amplitude Tl{t) at all times. 

Here we shall prove that, contrary to the scheme with Fock states, the scheme ( p7[p8| ) has an infinite statistical 
uncertainty for any time i > 0, in a sense to be precised below. As in the section with Fock states we introduce 

A{t) = (Barg : 0(i)|Barg : 0(t)) = |no|' eH^^H'. (50) 

As initial state at i = we take for simplicity a pure Bargmann state |Barg : 0(0)) with in average ||(?!)(0)||^ — N 
particles. Hq is taken as exp(— 7V/2), so that the initial state is normalized to unity. We now show that the mean 
value of A{t) is infinite for any t > 0. 

From (|^) and (|4|), it follows that the evolution of X{t) = ||0(<)||^ is described by 

dX = {dB\dB) + {dB\(j}) + {(j}\dB) = YKv^dt + d?7 (51) 

where 

dri = {dB\(l)) + {(l)\dB) (52) 
is a zero-mean noise term and its variance is equal to 



d?72 ^ ((dB|0) + {(}>\dB)Y = (0|dB)2 + {dB\4>)^ + 2{dB\(t)){(f)\dB). (53) 
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Inserting in this expression the noise correlation function (]48|), the first two terms on the right-hand side of (|52 
cancel. For the remaining term, a lower bound can be obtained by using the explicit expression ( ^9| ) for the noise and 
applying the Schwartz inequality ^: 

w = ^mmm = ^ ^ i^^r > dt, (54) 

r 

as in the previous sections, M is the number of lattice points. 

From (^), a lower bound on the increase of the moments of ( can be obtained: 



(55) 



/XV h 

and therefore on the moments themselves: 

(X")(i) > (X")(0) e^(")^ = A^" e''(")^ (56) 

Here we have set r = got/AVh and 

7(n) = n-|-n(n- 1)/M. (57) 
This leads to the lower bound for the mean value of A{t): 

(A)(t) - iHol' (expX) > iHol' ^-iV"e^(")^ (58) 

n=0 

As the ratio of two successive elements of the series tends to infinity in the limit n ^ oo for any < > 0, this results in 
an infinite value of (A)(t) and therefore of the variance of the statistical error on the many-body density operator. 

The same result could have been obtained by explicitly calculating the probability distribution for X at a given 
time t: as the evolution of the random variable Y = ln(X) is a simple diffusive motion in presence of a constant drift 
term, its probability distribution P(Y, t) has a Gaussian shape at all times t. By going back to the original variable 
X by means of P{X, t) = ^P{X^ the distribution of X is found to be equal to 

1/2 

X y^-KT 

It follows that the average of the exponential cxp(X) does not converge Sit X — +oo and therefore that the average 
of (|5^) is indeed infinite. 

One might wonder if the same conclusion holds for other indicators of the magnitude of the statistical error. For 
example, one could consider the mean value of the distance 5 between the Monte Carlo ansatz |Barg : 0) and the 
exact many-body state vector 

<5 = |||Barg:^)-|V')||^ (60) 
where ||. . .||^ is the usual norm in the many-body Hilbert space. As 

|||Barg:0)||^< |||Barg:^)-|V.)||^ + |||V.)|l^ (61) 
the mean distance {S) is bounded from below by 

{5)>(A''^)~\\m^. (62) 



Pi^X,t) = -(—\ e-|lNx/^)~(i-i,)rf. (59) 



^ For a fixed value of the norm of the field the sum J^^. |(/){r)|* is minimum if the density |(/)(r)p is a constant, which implies 

^^^V\4>{r)\*>Ut /MAV. 
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By a straightforward generalization of (|5^), we find that the mean value of A-"^/^ is infinite for any t > 0, so that the 
mean value of S is infinite. From the explicit form of the probability distribution (^9|) one can also check that the 
same result holds for any moment {S") of the distance 5, with an arbitrary real exponent a > 0. 

These results nicely demonstrate the difference between the convergence of the stochastic process for the field 
amplitude (j) and for the many-body ansatz. Since the deterministic term (^7|) is norm-conserving and the noise 
term (^) does not grow faster than one can apply a regularity theorem which guarantees that no finite time 
divergences can occur in the random dynamics of the field amplitude (f> and that all the moments of the stochastic 
wave function are finite (cf. chap. 4.3] and chap.6-7]). However, the statistical uncertainty is still divergent 
because the Bargmann ansatz involves an exponential function of the field amplitude ^(r) and not only a polynomial 
one as was the case with Fock states Hi. 



The simple scheme with coherent states 

In the previous subsection we showed that the statistical uncertainty for a simulation performed according to the 
Bargmann scheme is infinite. In we have identified among all the schemes based on coherent states the one (the 
so-called simple scheme) which minimizes the growth rate of the variance of the statistical error on the many-body 
density operator. Here we shall show that for any finite time the statistical uncertainty is infinite also for the simple 
scheme, similarly to the case of the Bargmann scheme. 

The scheme denoted in Q as the simple scheme with coherent states corresponds to a slightly more general ansatz 
than @: 

|coh : (0,n)) = He^^^^'^('")*'('")|0) (63) 

where the amplitude H(i) is now a stochastically evolving dynamical variable. The field amplitude cf) now solves the 
Ito differential stochastic equation 



dt 



d0(r) = ^ ho + gomr)f 0(r) + rfB(r) (64) 



with exactly the same noise term dB as in the Bargmann scheme, see (^9|). The evolution of the amplitude H is given 
by 

dn = -Il^AV(l)*(r)dB{r) = -Il{(t>\dB). (65) 

r 

The deterministic term in dcf) includes the usual mean field contribution. As it remains norm conserving, the results 
previously proved for the distribution and the moments of X{t) = ||0|p for the Bargmann scheme still hold for the 
present scheme. In particular, the stochastic differential equation (^) has a well defined solution according to the 
regularity theorem of pO|, H. Setting H = expS", equation (|65|) can be turned into 



dS = -{(f>\dB) -^{(P\dB)^. (66) 
Using the explicit expression for the correlation function of dB, this expression simplifies to 

dS = -mB) - ^ ^ AV^ |0(r)|^ . (67) 

r 

This equation for S can be integrated by quadrature proving the existence of a regular solution to (^5|) . 
As in the previous section we quantify the statistical spread in the simulation with 

A{t) = (coh : ((?!),H)|coh : ((/>,n)). (68) 

As the initial state at t = 0, we take for simplicity a pure coherent state |coh : (0(0), H(0))) with on average 
= ^ particles and an amplitude H(0) — exp(— Af/2), so that the state vector is normalized to unity. 



The time evolution of A(f) is given by 



dA=\\dB\?A = -^XAdt. (69) 
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Because of the nonlinear term proportional to X{t)A(t) in (|69|), the coupled system of ( p^ ) and ( |5l| ) does not fulfill 
all the hypothesis of the regularity theorem chap. 4.3] and therefore we can expect singular behaviors. 

In terms of the rescaled time r = got/AV H, the stochastic equation of motion for X{t) and A(r) can be rewritten 

as 



dX = XdT + dr] 
dA = XAdr, 

where the zero-mean noise drj satisfies, as in the previous section; 

d^ > —X^ dr. 

It follows from that 

A{t) = exp dr' X{t') 
and therefore its mean value can be written as the series 

oo 

(A)(r) = l + 5]A.(r) 

n=l 

whose coefficients D^ir) are given by integrals of the form 

Dnir) = [dTn dTn-1 . . . dT2 dn (X(t„)X(t„_i) . . . X(t2)X(ti)). 
Jr>r„>T„_i>...>r2>Ti>0 

Using the regression theorem of stochastic analysis |l^ and the fact that all moments ( p6| ) of X are increasing functions 
of r, it is easy to see that 

(X'=(T)X(r„_fe)X(T„_fe_i) . . . X(r2)X(ri)) (76) 
is itself an increasing function of t for any 1 < k < n and for t > T„_fe > Tn-k-i > • • • > > 0. This implies that 

(X(T„)X(r„_i)X(r„_2)...X(Ti)) > 

> {X{Tn-if X(t„_2) . . . X{ri)) > {X{Tn-2f ■ ■ ■ X {n)) >...> (X(ti)") > TV" e^(")^i . (77) 
For r > 0, we then have 



(70) 
(71) 

(72) 
(73) 
(74) 

(75) 



Dn{T) > 



(n-l)!7(n)" X 



e""" dw. 



(78) 



The integral in w can be explicitly performed and for values of n sufficiently large that ^{n)T > 1, it turns out to be 
bounded from below as follows: 



7(n)r 

e""' w"~^dw = {n-iy. 
which means that 



1 -e 



-7(n) 



n-1 



(7(n)T) 



>{n- 1)! [l - n(7(n)T)"-i e^'^^")" 



TV" g7(n)T 
7(71)" 



1 — n (7(»t.)t) 



Since 



lim — = +00 



(79) 



(80) 



(81) 



n— ►00 Efi 

the infinite sum of the En is infinite as well as the sum of the Dn which appears in (0). This proves that the variance 
of the statistical error on the many-body density matrix is indeed infinite also in the case of the simple scheme with 
coherent states. 

As for the Bargmann scheme, this result can be generalized to other indicators of the statistical error, e.g. the 
distance 6 between the Monte Carlo ansatz and the exact state vector 



<5=|||coh:(</>,n))- 1^)11^ 



(82) 



has an infinite mean for t > 0. 
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CONCLUSION 



We have shown that the simple Fock scheme can be rederived in a very elementary way using the Schrodinger 
equation for the state vector rather than the LiouviUe equation for the many-body density operator. Such a derivation 
is inspired from the derivation for fermions used in |8| . We have also given a proof of completeness of the Hartree-Fock 
state ansatz more detailed than in [Q. 

For the two schemes with coherent states considered in this paper, we have shown that the absolute value of the 
statistical error on the many-body state vector has an infinite mean for any finite evolution time. This negative result 
has been obtained by analytical means for a continuous time evolution. Of course, any actual numerical calculation 
is performed with a finite time step dt so that any given time t is reached with a finite number of steps. As the noise 
( ^ ) is bounded from above, this imposes an upper limit to the maximum value of the statistical error which can be 
achieved by any realization and therefore its mean value remains finite in the numerics. This is why this problem of 
the simple coherent scheme was not discovered in the numerical simulations of Q . 
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